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Abstract 

(N 

^ We consider supersymmetric models that include particles beyond the Minimal 

Supersymmetric Standard Model (MSSM) with masses in the TeV range, and that 
^ couple significantly to the MSSM Higgs sector. We perform a model-independent 

^ analysis of the spectrum and couplings of the MSSM Higgs fields, based on an 

effective theory of the MSSM degrees of freedom. The tree-level mass of the hghtest 
ON CP-even state can easily be above the LEP bound of 114 GeV, thus allowing for 

^ a relatively light spectrum of superpartners, restricted only by direct searches. 

The Higgs spectrum and couplings can be significantly modified compared to the 
^ MSSM ones, often allowing for interesting new decay modes. We also observe that 

^ the gluon fusion production cross section of the SM-hke Higgs can be enhanced 

with respect to both the Standard Model and the MSSM. 



1 Introduction 



Supersymmetry (SUSY) offers an elegant solution to the hierarchy problem, can explain 
electroweak symmetry breaking (EWSB) dynamically through renormalization group running 
and, if established experimentally, could open a window into the physics associated with length 
scales much shorter than can be probed directly. 

The most studied SUSY extension is the minimal supersymmetric standard model (MSSM). 
In the Higgs sector, the model incorporates two Higgs doublets, Hu and H^. After EWSB, 5 
physical scalars remain in the spectrum pQ. Assuming no CP violation, these can be classified 
as: two neutral CP-even scalars {h^ the lightest, the heaviest), one neutral CP-odd scalar 
{A^), and a charged Higgs pair {H^). The phenomenology in the Higgs sector is largely 
determined by the masses of these particles, by a mixing angle a that governs the relation 
between gauge and mass eigenstates in the CP-even sector, and by the ratio of the Higgs 
VEV's, tan/5 = Vu/vd- In the MSSM only two of these parameters are independent, and are 
conventionally chosen as niA and tan f3. Dependence on other sectors of the theory enters 
through radiative corrections, most notably in the mass of the lightest CP-even scalar, h^. 
This scalar is found to be below about 130 GeV [2], and together with the direct bounds 
imposed by LEP suggests the presence of some degree of fine-tuning in the MSSM. 

This has been one motivation to study extensions of the MSSM that can either relax the 
upper bound on the lightest Higgs mass [21 H], or alter its properties in a way that weakens 
the LEP bounds [5] . Often the extensions considered aim at addressing theoretical issues such 
as the /i-problem, i.e. how to link the supersymmetric Higgs mass parameter to the scale of 
EWSB, which in turn is set by the scale of SUSY breaking in the observable sector. Further 
theoretical constraints are also often imposed, such as the requirement of perturbativity up 
to a very high scale, the preservation of gauge coupling unification, or various simplifying 
assumptions that allow to more easily constrain the low-energy parameters. However, when 
exploring the Higgs collider phenomenology, it may be healthy to keep an open mind regarding 
such theoretical assumptions. 

If the constraints on the MSSM are to be relaxed, it is necessary to introduce new degrees 
of freedom (that interact with the MSSM Higgs sector) at or near the weak scale. It has 
been observed that even if the new particles in the Higgs sector are slightly heavier than 
the EW scale, the lightest Higgs boson mass may receive important contributions that can 
relax the LEP constraints [B] (see also |2| for examples with TeV scale vector-like matter). 
The new degrees of freedom may or may not be directly accessible at the LHC, but in either 
case their presence could potentially be inferred by studying the spectrum and couplings of 
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the lighter states. In particular, the observation of a SM-like Higgs with a mass significantly 
above 130 GeV, together with the observation of other superpartners, would provide a clear 
hint that the Higgs sector is more complicated than in the MSSM. 

We concentrate on supersymmetric scenarios with particles beyond those in the MSSM, 
under the assumption that they have order one couplings to the MSSM Higgs sector, and 
that they are heavier than, but close to, the weak scale. This allows to perform a model- 
independent analysis of the properties of the lighter states (i.e. those of the MSSM) by 
encoding the effects of the heavy physics via higher-dimension operators. As pointed out 
in [6] at leading order in 1/M (where M is the scale of the heavy physics) only two new 
parameters are introduced (corresponding to two operators in the superpotential). Therefore, 
even if the MSSM extension turned out to include a large number of degrees of freedom, their 
low-energy effects admit a rather simple parametrization. Furthermore, it was found that 
the 1/M effects can give rather important contributions to the mass of the lightest Higgs 
state. That the effects of such 1/M-suppressed operators can be as important as those of 
the renormalizable terms can be understood by considering the structure of the Higgs quartic 
couplings. In terms of the general two-Higgs doublet model (2HDM) parametrization of [9] 
(see Eq. ([s]) below), it is well known that of the seven independent quartic couplings only 
■^1) -^2 5 A3 and A4 receive contributions in the MSSM, at tree-level. The leading order higher- 
dimension operators contribute to A5, Ae and A7, so that they lead to qualitatively new effects 
in the Higgs sector. This is the same underlying reason that loop effects in the MSSM, which 
turn on all possible quartic couplings, can give sizable effects. Thus, the presence of the heavy 
physics allows a Higgs spectrum consistent with the LEP bounds, even if the SUSY breaking 
terms (in the top-stop sector) are of order a couple hundred GeV, thus alleviating the tensions 
found within the MSSM. 

Working still at leading order in 1/M, one finds not only contributions to the quartic Higgs 
couplings, but also higher dimension operators in the Higgs potential. These are essential in 
bounding the scalar potential from below. It was observed in [10], that taking these operators 
into account leads to the existence of new vacua that can be studied within the above effective 
field theory (EFT), and that these vacua have distinct properties. For instance, electroweak 
symmetry breaking is not necessarilly controlled by super symmetry breaking but, unlike in 
the MSSM, can occur already in the supersymmetric limit (this possibility was considered 
in the early days of supersymmetry [llj, though not in light of the EFT approach. See also 

"'^The earlier Ref. [5] also considered in detail the effects of higher-dimension operators on the MSSM Higgs 
sector in the context of low-scale supersymmetry breaking mediation. They also considered effects similar to 
those we include below, although they restrict to a study of the renormalizable part of the scalar potential. 
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Appendix B of Ref. [12], and more recently [I3]). These were dubbed "sEWSB vacua" (or 
supersymmetric EWSB vacua) in [10]. Typical features of the Higgs physics in the sEWSB 
vacua are order one tan /?, a heavy CP-even Higgs (H^) with SM-like properties, and relatively 
light charginos and neutralinos. 

In this work, we further consider the next order in the 1/M expansion. On the one hand 
these lead to the "first order corrections" to the physics in sEWSB vacua [lOj. Second, even 
in MSSM-like vacua their effects can be phenomenologically relevant. This can again be 
understood by referring to the quartic couplings Ai,2,3,4, whose size is set at leading order 
by the electroweak (EW) gauge couplings squared, hence are numerically small (the source 
of the lightness of the SM-like Higgs within the MSSM). The 0{1/M) operators from the 
superpotential do not contribute to these couplings, and therefore the 0(1/M^) effects become 
the leading order contributions from the heavy physics to the corresponding quartic operators 
in the scalar potential. Furthermore, up to numerical factors, the heavy physics gives a 
contribution of order f^/M^ which is comparable to the MSSM one for M ~ v/g, with g 
an EW gauge coupling and v a Lagrangian parameter of order the EW scale]^ At second 
order in the 1/M expansion several new operators appear. Within a given UV completion, 
the coefficients of these operators may or may not be related to the coefficients of the leading 
order operators (depending on the complexity of the MSSM extension). Here we do not 
impose any correlations, but vary the coefficients of the higher-dimension operators (in the 
superpotential and Kahler potential) independently. Our purpose is to survey the possible 
signatures in the Higgs sector, which could suggest the presence of an extended sector that 
might be more difficult to observe directly (for instance if it consists of SM singlets). 

Having developed the relevant formalism, we start in this paper a study of the Higgs 
collider phenomenology. We contrast the observed features against both the SM and MSSM. 
For instance, we observe a general enhancement of the gluon fusion production cross section 
of the SM-like Higgs, which is interesting at hadron colliders. Also noteworthy is the fact 
that the Higgs spectrum can be altered sufficiently to allow for new decay modes with rather 
significant branching fractions. Here we comment only on some of the possible signals, and 
defer a more complete study of the Higgs collider phenomenology to [H]. 

Higher-dimension operators in the SUSY context were considered in [T5l E] , and a more 
complete classification was presented in [TH], where field redefinitions were used to reduce 
the number of independent operators. The issue of the stability of the Higgs potential for 

^Note that sizable effects at order do not necessarily signal a breakdown of the 1/M expansion. In 

contrast, large effects at order in general signal such a breakdown. We take care to study parameter 

points where the next order effects are expected to be small by a simple criterion described in Subsection |4.2 
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MSSM-like minima (as opposed to sEWSB minima [10]) was considered in [T7], while the 
imphcations for fine-tuning in such scenarios was considered in [181 EH]- Higher-dimension 
operators can also have interesting consequences for the dark matter relic density [201 EI] as 
well as for cosmology ^22] and EW baryogenesis [231 IMl [2511261 [27] . 

This paper is organized as follows. In Section [2l we define the effective theory to be studied 
and work out the general expressions for the masses and couplings of the light Higgs degrees 
of freedom to order 1/M^. In Sectionjsjwe give simple analytic formulae in certain limits that 
allow to understand the qualitative features of the effective theory. In Section |4| we present the 
strategy to be used in the numerical study, to be undertaken in Section [5} There we comment 
on a selected number of features. A more complete study of the collider phenomenology 
will be presented in [Tl] . We conclude in Section |6} In Appendix |A| we consider several 
possible UV completions that illustrate how the higher-dimension operators in the EFT can 
arise. In Appendices |B| [C] and D we comment on custodial symmetry violation, give the 
chargino/neutralino mass matrices, and collect several Higgs trilinear couplings, respectively. 



2 Extended Supersymmetric Higgs Sectors 

We start by setting up the framework. Our point of departure is a generic supersymmetric 
theory with an extended Higgs sector, but where the degrees of freedom beyond those in the 
MSSM have masses of order M, assumed to be slightly larger than the EW scale. In this 
case, a model-independent effective field theory analysis is useful. We further assume that all 
supersymmetry breaking parameters are of order a couple hundred GeV, so that the heavy 
spectrum is approximately supersymmetric, with masses of order M. In this case, it is useful 
to write the effective theory in supersjTumetric notation, keeping track of supersymmetry 
breaking effects via a spurion superfield that gets a VEV in its F component (we do not 
consider D-term breaking here). In this section, we define the effective theory to be studied 
-which describes explicitly only the MSSM Higgs degrees of freedom- and work out the Higgs 
spectrum and their couplings. 

2.1 Generalized SUSY Two-Higgs Doublet Model 

The effects of heavy particles on the physical properties of the MSSM Higgs fields can be 
described by a tower of higher-dimension operators suppressed by powers of M . Our ultimate 
goal is to study the associated collider phenomenology, and we start by working out the 
spectrum and couplings of the light states (light compared to M). For the reasons spelled 
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out in the introduction, we work up to next-to- leading-order in the 1/M expansion. Next- 
to-next-to-leading order contributions are expected to be small, provided the 1/M expansion 



converges, a point we address in Subsection 4.2 



At leading order in 1/M, the superpotential reads 

W = fxH^H, + ^{HuH,f , (1) 

where HuH^ = H^H^ — H^H^ and ui is a dimensionless parameter that we assume to be 
of order one. Soft supersymmetry breaking can be parametrized via a spurion superfield 
X = msO"^, where sets the scale of SUSY breaking (note that we choose X to be dimen- 
sionless). Each operator in the superpotential (or Kahler potential) leads to an associated 
SUSY breaking operator through multiplication by the spurion, X. We will assume here that 
the coefficients of the SUSY breaking operators are proportional to those of the corresponding 
supersymmetric terms. Besides the Bn term, at leading order in 1/M one has 

M^spurion = ai^X(iJ„i7,)2 (2) 

for a dimensionless parameter ai, also taken to be of order one. This term leads directly to 
a quartic interaction among the Higgs scalar fields. 

At order 1/M^ there are no operators in the superpotential, but several operators enter 
through the Kahler potential: 

K = Hle^^Ha + Hle^^Hu + AK^ + AK'', (3) 



where 



M 



iMf ' W " 



(i/„i7,)+h.c. (5) 



We separated the higher- dimension contributions in the Kahler potential into those that 



violate the custodial symmetry, i^!^, and those that respect it, K'^ (see Eqs. (32) and (33) 
below, as well as Appendix |B|). 
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The above dimension-6 operators lead to associated SUSY breaking operators as follows: 



i 



Cl 



spunon 



|M|2 



[(7iX + 7i*Xt)+/5iXtx] {Hle^'^E, 



+ 



+ 



_C2_ 

m 

C3 



\M\ 



[(73X + 73*Xt) + /33XtX] {Hi e'''H^){H\ e^^if, 



(6) 



spunon 



C4 

IMP 



+ 



[(74X + 7:xt)+/?4XtX] 



C6 



IMP 



■(76X + ^eXt + /JeXtX)//] e^^iJ, 



(7) 



+ nn^(77X + + /57XtX)i7t g^v^i/^ 



IM 



{H^Hd) + h.c. 



Here we used the fact that operators like XH\ e'^^Ha or XH^ e^^Hu can be set to zero by a su- 
perfield redefinition [16]. We did not write explicitly the usual soft breaking masses m|^^ and 
m|^^, which correspond to operators of the form X^^XH^ e'^^Hu and X"'"Xif^ e'^^Hfi, although 
such terms are understood. We assumed that each type of soft breaking operator vanishes in 
the absence of the corresponding SUSY preserving one, i.e. we define the parameters 7^, 6i 
by factoring out the associated q. This will be the case if the corresponding statement holds 
in the UV theory that induces the higher-dimension operators we are considering, and is a 
property of several realistic SUSY breaking mediation mechanisms (in simple extensions one 
obtains a strict proportionality, as illustrated in the examples described in Appendix |A]) . No- 
tice also that operators like X^HuHdY or X'^ X^HuHdY in the Kahler potential are equivalent 
to the superpotential operators of Eqs. ([T]) and ([2]), and since we are taking uji and ai as free 
parameters, there is no loss of generality in omitting them from Eq. ([T]). 

It is straightforward to work out the scalar potential that follows from Eqs. ([T|, ([2]), ([3]), 
^ and ([7]) which, at the renormalizable level, takes the form 

= mlHlH^ + m^H^Hd - [bH^Ha + h.c] 



(8) 



where we added the soft masses , 



|r and b, and m?. 



u,d 



m|^^^ + l/ip also include 
the supersymmetric mass term, |/ip. The quartic operators are defined in such a way that 
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the Xi,i = 1, . . . , 7, map exactly to those used in [9] with their alternate convention for the 
hypercharges of the two Higgs-doublet fields. In the MSSM limit, the quartic couplings are 

,2 



(0) 



A 



(0) 



A 



(0) 



r 

' 2 



(9) 



with the rest vanishing. At order 1/M (henceforth referred to as dimension-5), only A5, Xq 
and A7 receive contributions: 



AA 



(5) 



M 



AA 



(5) 



AA^' 



(5) 



UJi 



M 



(10) 



At order all the quartic couplings receive contributions: 



(6 



AAi' 



AA 



(6 



(6 



AA 



AA^ 
AAf 



■2(c, + c.)A_+4cA-^-2c,ft-^, 

2 

n 777, II 

-(Ci + C2 + 2C4)|^ + 2(06^6 + Ct'^t)^ 

2 

-(ci + C2 + ^cs)— + 2{cQde + CjSrjj^ 



M2 



2(C676 + C777) 



M2 



■(^6 + 2c7)^ + (2ci7i + C373 + C474)^ - CePejfi 

■(^7 + 2C6)^ + (2C272 + C373 + C474)^ " ^7/^7^ 



where we assumed, as we will do for simplicity in the rest of the paper, that all parameters are 
real. At this order we should consider also dimension-6 operators in the scalar potential that 
give parametrically comparable effects. The dimension-6 operators generated by integrating 
out the F- and D-terms take the form 



non— ren. 



1 

M2 



H^H,\%XsHlH, + X'.HlH^))] + Xg\H^H,\^ + X,o{HlH, 



+ XniHiH^f + X,2{HlH,){HiH^)) [H^H, + HIH^] (12) 
+ Xi:,{HlHdf + Xi^{HlHdf{HlH^) + Xi,,{HlHa){HiHuf + Xi,{HlHuf \ , 
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where 



As =UjI- (Ci + C3)y , Ag =UjI- (c2 + Cg)^ , 

Ag = -(C6 + Cy)— , Alo = C6 , 

An^ci^, A.Hca + c.)(^, (13) 

, (^7^ + g") , g' {g'-g") 

Al3 = Ci ^ , Ai4 = C3— - Ci- 



4 ' ^^2 ^4 

, g' jg' - ^7^^) , jg' + ^7^^) 

Al5 = Csy - C2 , Ai6 = C2 . 

The operators with coefficients Ag and Ag are essential in stabihzing the sEWSB vacua dis- 
cussed in [TU]. At order higher-dimension operators involving two derivatives need also 
be included, since after EWSB they lead to contributions to the Higgs kinetic terms. After 
canonical normalization, these give additional contributions to the masses and couplings of 
the Higgs fields. These operators are found to be 



C D 



M2 

+ C2 {[{D''Hu)^Hu] {HlHu) + [{D^Hu)^Hu][{D''Hu)^Hu]) 

+ C3 (^[{D^Hu)^H^] (HlHd) + [{D^HdYHd] {HIH^) + 2[{D^Hu)^ Hu][{D^Hd)^ H^] 
-c^d^{HlHl)d^{HM + {[c,{D'H,yH, + cj{D^H^yH^] + h.c.) | ,(14) 



where = D^D^^ and = df^ + igW^ + ig'YBfj, is the gauge covariant derivative {Y = +1/2 
for Hu and y = -1/2 for Ha). 

Finally, although there can be 0{1/M'^) contributions to the gauge kinetic terms, we do 
not consider them here (but see [8]). 

2.2 Higgs Spectrum 

We obtain the spectrum and couplings of the Higgs fields in two steps: ffist we neglect the 



non-canonical kinetic terms that follow from Eq. (14), minimize the potential (which does 
not require canonical normalization), and obtain general expressions for the (at this point 
unphysical) masses and couplings in terms of 

A, = Af ) + AAf ^ + AXf + AAl-^°°P , (15) 
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where the last term represents the 1-loop corrections (to be included in the numerical analysis 



of Section 5|, which will be discussed in Subsection 4.5 In this section we concentrate on 



the tree-level contributions [i.e. the first three terms in Eq (15)]. Second, we make a field 
redefinition to arrive at canonical normalization in the presence of the operators of Eq. (14), 
and rediagonalize to obtain the physical spectrum. The first step defines CP eigenstates 
h^, H^, A^, H^, with mass parameters m^o , ""^^o , ""^^o , "^^i as follows, 

( G-o 



+ 



? ? , (16) 



(17) 



where sp = sin/3 and C/3 = cos/? with tf^ = Vu/vd and v'^ = vl + v"^ ^ (246 GeV)^. For small 
fluctuations, G^'^ are the eaten Goldstone bosons that will be set to zero in the following. 
Minimization of the potential requires 

■Sfl [3Ai6 + 5Aiito^ + 2Agta^ + SXgtn^ + Agt^^ + Aiot«^] , (18) 



4M2'^ 

ml = 6t/3 - yt/3S^ [A7 + Ast^^ + SAet^^ + Alt 

^ ^ ^ 

- ^tp4 [^11 + ^8^^' + SAgt^' + 2Xst/ + 5Aiot^^ + SX^t/] , (19) 



where A3 = A3 + A4, A3 = A3 + A5, Ag = Ag + A14, Ag = Ag + A15 and Ag = Ag + A12. The masses 
of the CP-odd and charged Higgses, and if^, can then be written as 



f\ 0\ -1 \ -2\ "^^ 

AO = ^ - -jtp [X7 + 2X^tp + XQtp ) - ^2 



2 

m — 



tp (A7 + 2A5t^' + Agt^') - tpsl [An + Agt^^ ^ ^^^^-4| ^ ^20) 

2 nA 



= m\o + - (As - ^4) - ^ si (A^ + 2X^rp^ + Agt^^) . (21) 

The mass matrix for the CP-even states in the (/i^, hy) = {caH^ — SahP, Cah^ + SaH^) basis is 

1 ~^i3^ j _|_ 2 2 I "^^ ~'~ ^-^6^/3^ + Ait^^ A7 + Ast^^ + Xet 
^13^ ^13^ j V A7 + Ast^^ + Agt^^ A2 + 2A7t^"'^ + Asl^ 

f^s^ / Agt^"*^ + Agt^^ + 4Aiot^^ + SAist^^ An + Agt^^ + 2Agt^^ + Agt^^ + Aiot^ 



+ 



\^ An + A'gt/ + 2Agt^2 ^ ^^^-3 ^ ^^^^-4 3^^^ ^ ^^^^^^i ^ ^,^^-2 ^ 



which leads to mass eigenvalues and a mixing angle a given by 



m 



HO,hO 
C2a 



TtM'±J (Tr MY - 4 det 



12 



{TiM^y -AdetM'^ 

2 
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Ml -MI2 



{TiM^f -AdetM^ 



(23) 
(24) 

(25) 



In order to find the physical masses at order we need to include the operators in 

Eq. (14), that contribute to the kinetic terms of H^, A° and H^. Canonical normalization 



is achieved by the field redefinitions 



where 
A, 



El 



2 
1 

2^ 



- (1 - a' , 



(26) 



2M2 



2M2 



2M2 



2M2 



■^^ 1*^2(1 + C2a) + Citp^il — C2a) + -^C^ [s^^ — (l — t^^) C2a — 2t^"'^S2a] 

+ C7 [-S2a + 2t-p\l + C2„)] + C6 t^' [2(1 - C2„) - ^^20] ] , 

4 I (C2 - Clt;^^) S2a - ^C3 [(l - t^^) S2a - 2t^^C2a\ 

+ CiS^'^Ca+fjSc+p + C7 (C2a + 2t^^S2a) + Ce (-2S2a + t^^ C2a) } , 

4 1^2(1 - C2„) + Cit^2(l + C2a) + ^C3 [s^' + (l " t/) C2„ + 2t^ls2j (27) 

+ C4S^^S^+^ + C7 [S2a + 2t^^(l - C2a)] + Cq [2(1 + C2a) + t~^^ S2a\ ] , 

I ^(Ci + C2) [1 - 4 (1 - 14t^2 + t-p% + ^C3 [3 + 4 (1 + 2tf + t^^)] 

+ C4 + 2t^i4 [c6 (1 + 4) + cr (1 + tfsl)] } , 



(3s/ + 1 - + t/) + t/(c6 + Cyt/) + ^(Ci + C2)tJ 
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At this point, one needs to perform a further rotation (by an angle 7) in the CP even sector, 
that rediagonahzes the corresponding mass matrix: 



where, to order 



t 



27 



(28) 



m|^o(l 



We do not expand the denominator to cover cases where m^o 
squared masses are finally given by 



(29) 

The physical CP-even 



mlo (1 - Ai) + m^o (1 - Di) + 5i(m^o + m^o)c^s^ , 
niffo (1 - -Di) ci + mlo (1 - Ai) - i?i(m^o + mlo)c^s^ , 



(30) 



while the physical CP-odd and charged Higgs masses are given by m^o(l — -E'l) and m^±(l — 
Fi), respectively. Whenever 7 ^ 1 we have ^ m^o (1 ~ ^1) and ^ m^o (1 — Di). 



However, there are regions where the denominator in Eq. (29) is small and all orders in 7 
should be kept, in spite of it formally being of order 

There is one additional source of corrections at order that affect the spectrum, 

as well as the Higgs couplings to be discussed in the next subsection. The two-derivative 



operators in Eq. ( 14 ) give contributions to the gauge boson masses as follows: 

,2 



ml 



4 



2 2 
9zV 



m 



w 



m^c^ (1 + aT 



where 



aT 



2M2 



(31) 
(32) 

(33) 



is the contribution from the higher-dimension operators to the Peskin-Takeuchi T-parameter. 
Note that, as mentioned in Section 2.1, only the operators in Eq. (|4]) contribute to T. At loop 
level, there are other contributions to the T-parameter from the distorted Higgs spectrum 
(heavier SM-like Higgs in addition to mass splittings among the non-standard Higgses), as 
well as from custodially violating mass splittings in the superpartner spectrum [32]. For the 
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time being, we note that keeping mz fixed at the observed value, Eq. (31) imphes a shift in 
V at order 1/M^: 



At; 

V 



2M2 



(34) 



so that the VEV v in this model differs from the SM value of 246 GeV by terms of order 
l/M^. At this order, only the terms in Eq. (21) and (22) that are proportional to Ai, A2, A3 



and A4 -which are non- vanishing in the MSSM limit- are affected by the shift in the VEV (in 
particular, m\ is not affected). In practice, at this order in the 1/M expansion it is sufficient 
to use w = (1 + /\v/v) X 246 GeV in Eqs. ([l6])-([30]), with /\v/v given in Eq. ([34]). 



2.3 Higgs Couplings to Fermions and Gauge Bosons 

Here we focus on the couplings of the Higgs scalars to gauge bosons and fermions. Recall 
that in the MSSM the couplings of the neutral Higgs bosons to fermion pairs, relative to the 
SM value, grrif /2mw, read (for up-type and down-type quarks) 

hit : Ca/s/3 Hit : Sa/s/3 Ait : l/t/3 

- - (35) 

hbb : —Sa/Cj3 Hbb : Ca/cp Abb : t/^ , 

while the couplings to a pair of gauge bosons relative to the SM value {grriw for V = W and 
gmz/cw for V = Z) are 

hVV : Sf3-a HVV : C/3-„ . (36) 
These couplings are changed by the higher- dimension operators in two ways: indirectly 



through the mixing angle a [see Eqs. (22), (24) and (25)], and directly due to the effects 



associated with the operators of Eq. (14), which include the field redefinitions of Eqs. (26) 



The latter involve a rescaling and therefore cannot be parametrized as a rotation by an ef- 
fective angle. These effects correspond to the mixing of the MSSM Higgs fields with heavy 
degrees of freedom that are not included explicitly in the effective theory, and appear first at 
order l/M^. 

Since it is convenient to give the Higgs couplings relative to the SM values, we need 



to take into account the shift in the Higgs VEV, Eq. (34), which implies a shift in the 



Yukawa couplings relative to the SM values: Ayf/i/f = —Av/v. Together with the Higgs 
field redefinition of Eqs. (26) these induce a shift in the Higgs-// couplings at order l/M^ 



compared to Eqs. (35) (on top of the shifts implicit through the mixing angle a). 

Similarly, the normalization of the Higgs-gauge boson couplings to the SM value needs 



to take into account the shifts in the gauge boson masses given in Eqs. (31) and (32). In 
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addition, the two-derivative operators in Eq. ( 14 ) also give direct corrections to these vertices 
from the terms quadratic in the gauge fields and linear in /i° or H^: 



where 
Sghzz 

Sgnzz 

^ghww 



2M2 



2M2 



{4S^ (C2C« - Cit^^S„) + si [C7 {2Ca+f3 + Ca-/?) + Ce^^^ {2Ca+p - Ca-fs)] } 



5g 



HWW 



|2S/3 [C2Ca - Clip Sa) + 2C3)f:^ SpCa+f} 
+ si [C7 (2Ca+/3 + Ca-p) + Cet^^ (2Cq+/3 - Ca-p)] } , 

{2s^ (c2Sq, + Cit^^C„) + 2c3ti^^slSa+f3 
+ sl [Cj {2Sa+f3 + Sa-is) + Ce^/ (2s„+^ - Sa-p)] } ■ 



(37) 



These should be combined with the field redefinitions Eqs. (26). Note that the corrections to 
the Z and W couplings are different only due to the custodially violating coefficients Ci, C2 
and C3. 



Putting all these ingredients together, the couplings of Eqs. (35) and (36) are then gener- 
alized as follows. For the light CP-even Higgs, h^, we get 



htt 
hbb 



1 



-a+7 



1 - 



Av 1 



1 

C/3 



'a+7 



V 2 
Av 1 



Ai - -BiSa-^ + -{Di - Ai)scS^ 



V 2 
Av 1 



1 -Ai + -BiCa-^ - -{Di - v4i)c„s. 



hZZ : Sf3_a-'r ^ ~ ^ ~ ^i)'^/^-"'^^ 



hWW : s^_a-7 1 + 



V 2 
+ c^Sghzz - s^6gHzz 
Av 1 

~ ~ 2 

+ C^^ghWW — S-i^gHWW 1 



(3J 
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while for the heavy CP-even Higgs, H^, we get 



Htt 



Hbb 



HZZ 



HWW 



1 

1 

C/3- 



(■ 


Av 


V 


(■ 


Av 


V 


1 + 


Av 


V 



- -Ai 
2 



-BiCa-^ - -{Di - Ai)SaC^ 



:Ai 



1 



- -{Di - Ai)caC^ 



1 



{Di - Ai)Cf3-aC^ 



(39) 



+ c-y5gHzz + s-y5ghzz , 
At; 1 



C/3-Q-7 I 1 + 



V 2 



aT - -Ai 



Q+7 



-iDi - Ai)Cfj-aC^ 



where a is determined by Eqs. (24) and (25). The couphngs to fermion pairs of the CP-odd 

and charged Higgses, A^ and H^, which are independent of a, are obtained from those in the 

MSSM by multiphcation by (1 — Av/v — ^Ei) and (1 — Av/v — ^Fi), respectively. Here 7 

is defined by Eq. (29), Ai, Bi, Di, Fi and Ei are as defined in Eq. ([27]), Av/v is given in 

1^ 



Eq. (34), aT is given in Eq. (33j)|jand dgnvv and Sgnvv are given in Eqs. ( |37| ). The hVV and 
HVV couplings in Eqs. (38) and ([39]) are given relative to gniz/cw and gm\Y respectively, 
where mz and my/ include the order corrections as in Eqs. (31) and (32). The Higgs 

couplings to fermion pairs are given relative to grrif /2mw, which is the SM Yukawa coupling. 

Trilinear interactions involving one gauge boson and two Higgses are also of great phe- 
nomenological relevance, and receive corrections from the higher-dimension operators of Eq. ([3]): 



^ {{vzhaH^A^ - r^zAHA^'d^H') - {vzkAh'>d^A' - VzAhA'^d^h')} 
VzH+H-Z^ {H+d,H- - H-d,H+) 



2cw 
igc2w 



2cw 
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- {r]w±HHTH^d^H - -qwiHTHH d^H^) ] + h.c. 
{vw^HTAH-d.A"^ - vw^AH^A^d.H-} + h.c. 



(40) 



We give the detailed form of the coefficients rjzhA, VzAh, Vzha, Vzah and rjzH+H-, as well as 
those for the interactions with a single W, in Appendix ]D] 



Eqs. (38 1 and (39 1 are the tree-level expressions. There are loop level contributions to aT, discussed in 
Section |4.4[ that have to be added. Since experimentally aT is bound to be small, one can neglect it in hWW 
and HWW for the phenomenologically allowed points. 
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3 Generic Features 



In the previous section we presented the effective theory in the Higgs sector up to 0{1/M^), 
where M is the scale of the heavy physics. It is useful to obtain simple analytical expressions 
that hold at order 1/M, since these will determine the qualitative features induced by the 
heavy physics. In this section we perform such an analysis, and consider several limiting cases 
that clarify the generic features of the numerical study to be undertaken in the next section. 

We start by considering the masses of the CP-even Higgs states, m^o and m^o, as given 
in Eq. (23) with Eqs. ^ and (10). Formally expanding to leading order in 1/M one has 

\Dim-5 



[m 



where the well-known MSSM contributions, (m^ 

.2 )Dim-5 (Am^o 



corrections due to the heavy physics, (Am^o 

1 
2 

1 



,o;^^^^ and (m^o; 

2 ^Dim-5 



(41) 

MSSM^ and the leading order 



are given by 



>MSSM 



(42) 



{Ami ^^'"^'^ 



''hO,HO) 



-UJiV 



2s 



2/3 



ii 

M 



1 ± 



m, 



rrir 



D 



aims 
M 



It 



mr 



D 



provided one has the inequality 



2s 



2/3 



m) ^ 



m 



2/3 



aim^ 
M 



mr 



(43) 



D = 



m\ + miy -Acli^m\m% 



For the opposite inequality, the expressions for m^o and m^o in Eq. (42) should be inter- 
changed. For typical values of the parameters we are interested in, the inequality Eq. (43) is 
generally obtained for large while the opposite inequality results for small m^. Near the 



point where the two sides of (43) coincide, the formal expansion to order 1/M can be rather 
inaccurate due to the near degeneracy of the two states (the mixing angle a can receive large 
corrections). Nevertheless the explicit expressions in Eq. (42) are useful, for instance showing 
how the sign of the correction to the MSSM result depends on the signs of cui/i and uiai. 

It is interesting to consider the t/j dependence of the correction to m^o and m^o- In the 
following we assume that the inequality Eq. (43) holds. As discussed above, the complemen- 
tary region can be obtained by simply interchanging and in the relevant statements 
below. For of order one (S2/3 ~ 1, C2/3 ~ 0) one has 



(Amio 



2 \Dim-5 



-UJiV 



4yU — aims 
M 



(44) 
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while for large (i.e. S2/3 ~ 0, |c2/3| ~ 1) one has 

„ „. . I tha > mz , 

' \ -^it^' (^) < . ^ ' 

Thus, the corrections are most important in the small tp regime (if ^ rnz), and this region 
will allow to more easily evade the LEP bound on m/jO. This is a well-appreciated feature 
in models extended by singlet fields, such as the NMSSM and its relatives. We can easily 
estimate how large m/jo can be for ^ 1, where the MSSM contribution vanishes. Setting 
«! ~ — 1 and fj, ~ nis, we have 

rriho ~ X C(5/i/M)i/^ . (46) 

v2 

For ~ 1 this can easily be above the LEP bound. For the case of large tf3 the dimension-6 
corrections can play a crucial role in lifting the lightest Higgs mass from the MSSM limit, 
as will be illustrated in the numerical analysis of the next section (notice that dominant 
dimension-6 effects do not signal a breakdown of the EFT analysis if the dimension-5 contri- 
butions are small due to a t/3 suppression or an accidentally small coupling ui). 

It is also interesting to consider the decoupling limit. Expanding in powers of m|/m^, 
one gets 



i^rnlor-' - ^^iv^s,, (^ ^- J , (47) 

which shows that in a large region of parameter space the leading-order correction to m/jO 



is only suppressed by about S2/3 compared to Eq. (46) [in general, this correction has to be 
added in quadrature to the MSSM result, according to Eq. (41)]. 

In Fig. [1] we show the maximum tree-level value of the lightest Higgs mass ruh as a 
function of rriA, for small tan (3, fixed representative values of the dimensionful parameters, 
and assuming that all dimensionless parameters take values at most equal to one in absolute 
value. We show the tree-level value of rrih up to order 1/M, which is obtained for ui = —ai 
(dashed-dotted curve). We also show the maximal values of rrih up to 1/M^ effects (dashed 
curve), and the MSSM (solid) curve for comparison. We see that the effects of the higher- 
dimension operators can be quite substantial. We emphasize that the effects of order 
or higher are expected to be much smaller, as mentioned in the introduction and discussed in 
Subsection 14. 2[ 

For the heavy CP-even state we find that (Am^o)^'™'^ ~ —^Uiv'^aims/M for tan/5 ^ 1. 
For tan/5 > 1, (Am^o)^™"^ ^ -Uiv'^aims/M if > mz, and (Am^o)°™"^ ^ if < 
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Figure 1: Lightest CP-even Higgs tree-level mass, nih, as a function ofniA, for tan 13 = 2. The 
dashed blue line corresponds to the maximum value ofm^ to (9(1/M^) when the dimensionless 
coefficients of the higher- dimension operators are allowed to he as large as 1 (in absolute 
magnitude). The dashed-dotted red line corresponds to the maximum value of to 0{1/M) 
under the same assumption. The solid green line corresponds to the tree-level MSSM result. 

mz- In the decoupling limit m\ ^ m|, we get 

(Amlor-^ = -\co,v^[S + cosm] [^) . (48) 

Hence, the heavy CP-even state receives corrections at the leading order in the 1/M expansion 
only due to the SUSY breaking operator, Eq. 

Similarly, the charged Higgs mass takes the form 

2 2 I 2 2 (^I'^s 



T^H± = m^ + m^- -UiV [-j^J , (49) 

and, at leading order in 1/M, gets corrections only from the SUSY breaking operator, Eq. 

Lastly, we consider the corrections to the Higgs couplings at leading order in 1/M. We 
start from the mixing angle a which, at this order is given by 

^ = i ^A + ^z \ ^ uJiv^ {-'2fi{m\ - m|) + S2i3{m\ + m|)aim,) 
^" ^^\m^ — m|/ Mc2/3(m^ — m|)2 



In the decoupling limit this simplifies to 



t2a = ^2/3 + 2t2(3X + _2 „ 77 -2/i + S2/3aim^ , (51) 
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where x = m%/m\. Taking into account that a is in the fourth quadrant, while /3 belongs to 
the first one, one gets 

a = p - - + S2(3C2f3X + -2ij + S2f3aim, 

= /3 - ^ + A^^^x + A(2)/M . (52) 



This implies that the couplings of the CP-even Higgs fields to two gauge bosons [see Eqs. (36) 
are 

hVV : l + C(x2,fVM2) , HVV : A^^^x + A(2)/M . (53) 

Note that the couplings of the light state to gauge bosons do not receive corrections at order 
1/M and are not expected to deviate very much from the MSSM ones, while those of the 
heavy state are expected to get larger corrections. The couplings to fermion pairs relative to 



the SM are [see Eqs. (35) 



htt : l + t-i(A«x + A(2)/M) , Htt : -t;^^ [l - t^(A«x + A(2)/M)] , 

hhh : 1 -t^(A«x + A(2)/M) , Hhh : -t p [l + t-p\A'^^^ x + A^'^^ M)] . ^^^^ 

Thus, in the decoupling limit there could be important variations of the light state couplings 
to the up and down sectors with respect to the SM predictions. However, the hit coupling 
remains SM-like in the large tan/5 regime. Similarly important variations occur in the cou- 
plings of the heavy Higgs to the up-type and down-type fermions, except for the Hhh coupling 
in the large tan (3 regime, where the variations are small. 



4 Numerical analysis: Parameters and constraints 

In the previous sections we worked out the spectrum and couplings of the Higgs sector in a 
softly broken supersymmetric theory, under the assumption that there is a set of particles that 
couple to the MSSM Higgs fields, but that have masses parametrically larger than the weak 
scale. We also assumed that the SUSY breaking mass splittings in the heavy supermultiplets 
are small, so that their masses have a nearly supersymmetric origin. In this section we 
define the regions or parameter space and discuss several constraints that will be used in the 
numerical exploration of the effects of the heavy physics on the MSSM Higgs sector, performed 
in Section [sl and further expanded in [T3] . 
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4.1 Parameter space 

We start by defining our region of parameter space. We assume tliat all SUSY breaking mass 
parameters (as well as /i) are of order the EW scale (a couple hundred GeV). For definiteness, 
we take /i = = 200 GeV, where is the F-component of the spurion superfield. We 
also set M = 1 TeV. In addition, we scan over the dimensionless parameters defined in 
Eqs. 0-0 as follows: 

• \uji\, |ci|, |c2|, Icsl, |c4|, Icel, IctI e [0,1]. 

• lAI, |7.|, e [1/3,1] for z = 4, 6, 7. 

Recall that the coefficients uji, ci, C2, C3, C4, Cg and C7 set the size of the 1/M and 1/M^ 
suppressed operators (both SUSY-preserving and SUSY-breaking). To be definite we assume 
that their values are at most one in absolute magnitude, but it should be clear that if they 
were larger the physical effects would be correspondingly larger. Our choice for the ranges of 
|«2|, l^il refiects our assumption that the SUSY-breaking operators are proportional 

to the corresponding SUSY preserving ones, so that these are parameters of order one. 

We consider two representative values of tan/3: tan (3 = 2 and tan/? = 20. We also vary 
the CP-odd mass up to 400 GeV, which is still below M, ensuring a proper separation between 
the light and heavy scales, as required by the EFT analysis. 

Note that scaling /i, and M by a common factor leaves the corrections due to the higher- 



dimension operators, Eqs. (10) and (11), unchanged [though not those of Eqs. (12) and (14)]. 
In particular, the leading order (dim-5) operators depend on fi, rris and M only through the 
ratios fi/M and m^/M. Thus, these effects could be relevant even if the scale of new physics 
is much higher than we envision here, if SUSY breaking and/or fi are correspondingly larger. 
Even though we do not scan over the values of /i, and M (but rather fix them as specified 
above), our results should be qualitatively applicable when all these scales are higher, keeping 



the ratios fixed [the difference arises at 0{1/M ) through the operators of Eqs. (12) and (14)] 



4.2 Uncertainty from higher orders and the EFT expansion 

As was mentioned in the introduction (see also Fig. [T]), the contributions of order 1/M and 
1/M^ can be phenomenologically sizable: the 1/M effects turn on Higgs quartic couplings not 
present in the MSSM at tree-level, while the 1/M^ effects can easily be comparable to the 
MSSM contribution which is set by the weak gauge couplings. However, one should make sure 
that the next order can be reasonably expected to give a small contribution, since otherwise 
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it could signal the need to resum the 1/M effects to all orders (in which case the details of the 
UV completion are essential and the EFT approach ceases to be useful). It can also happen 
(and does happen in our random numerical scans) that there are accidental cancellations 
between the MSSM contributions proportional to and those at order [see Eq. (11)], 

or also, if uoi is somewhat suppressed, between the order 1/M and 1/M^ effects. In such 
cases, the next order corrections can have a larger impact than naively expected, and our 
1/M^ analysis would fail to capture the quantitative properties of such points in parameter 
space. 

We perform a simple test to assess the robustness of a given point in parameter space 
against higher-order corrections that have not been computed (and that would depend on 
new coefficients that are arbitrary from the EFT point of view). Since the most important 
effects are expected to enter through the Higgs quartic couplings, we model the next order 
corrections as follows]^ keeping all the Lagrangian parameters fixed (m^, m^, 6, cui, Cj, ai, 7^, 
67, for i = 1,2, 3,4, 6, 7), we modify the quartic couplings by hand as: 

Ai ^ A, ±2Max{|u;i|, |ci|, Icsl, Icsl, |c4|, IcqI, IC7I} (^^j , i = l,...7, (55) 

i.e. we allow for 1 /M'^-suppressed operators with dimensionless coefficients as large as those 
of the leading two orders in the 1/M expansion [the factor of 2 models numerical factors 
that may appear, as we have seen in the order 1/M^ expressions, Eqs. (11)]. We then solve 
the minimization equations, (18) and (19), with these new Aj's, which leads to values of v 
and tan (3 that are slightly different from their values in the absence of the modification. 
The amount by which these two observables change should give a reasonable measure of the 
sensitivity of the given parameter point (truncated at order 1/M^) to effects from the next 
order 1^ We restrict to parameter points for which the above procedure leads to a change of no 
more than 10% in v. This should be taken as no more than an order of magnitude estimate 
of the uncertainties associated with the truncation of the tower of higher-dimension operators 



''Notice that at order there are only two new operators: the superpotential term {HuHdY and the 

associated SUSY breaking operator. There are no new Kahler operators at this order. However, there are 
effects associated with the lowest order coefficients, e.g. proportional to wiQ. 
^Notice that the way to interpret a VEV different from 246 GeV in the presence of the modification 



Eq. ( |55| ) is as follows: one should rescale all mass parameters by an appropriate factor so that one recovers 
the "observed" value for v. This corresponds to normalizing to the measured Z mass, and does not change the 
physics since all mass ratios are kept fixed. Since all mass scales, in particular the physical Higgs masses, are 
rescaled by the same factor, we see that the change in the VEV due to the modification Eq. (551 corresponds 
directly to a change in the spectrum. There can be additional contributions to the spectrum from the higher- 
dimension operators, not associated with the shift in the VEV, but we expect that these are of the same 
order as the effect from the shift in the VEV we have described. We therefore consider the above a reasonable 
estimate of the uncertainties associated with the higher order operators. 
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at order We impose a looser constraint on tan/5: for points with tan/3 = 2, the point 



is allowed if 1.5 < tan/3 < 2.5 after the modification Eq. (55); for points with tan/? = 20, 
the point is allowed if 15 < tan /3 < 25 after the modification of the A's. The rationale for a 
looser constraint in the shift in tan (3 is that these shifted values would still be representative 
small and large tan/3 cases, respectively, and no dramatic change in the Higgs properties is 
expected within the above ranges. 

4.3 Global versus local minima 



The potential given by Eqs. ([8j) and (12), V = Vrcn. + Kion-ren., has in general several minima, 
which may also include CP-violating or charge breaking VEV's. Without loss of generality, we 
can parameterize all minima by {Hy) = (0, f„) and {Hd) = {vd + ivs, vcb), where v^, vs and 
VcB are real. We can choose this form for (Hu) by performing an appropriate SU{2)l rotation. 
It is also clear from the fact that the scalar potential depends only on H^Hu, H^^Hd, and HuHd 
that, having set if^ = 0, it depends only on = vcb^ so we can assume that this latter 
VEV is real and non-negative. We look numerically for all solutions of the minimizations 
conditions d^^V = d^^V = d^^V = d^^gV = 0, which are polynomial in the variables, in order 
to check that the minimum being considered is the global one, and in particular preserves 
charge and CP. In this work we do not consider the possibility of explicit nor spontaneous CP 
violation. Also, we do not consider cases where the minimum is metastable but long-lived, 
although such a case can in principle occur and could be phenomenologically viable. 

We also recall here the possibility of having different types of (global) minima, as em- 
phasized in [To]. The point is that at the level of quartic couplings, the Ag and A7 operators 
in Eq. ([8| lead to runaway directions. These are stabilized by dimension-6 operators in the 



scalar potential, Eq. (12). As a result there can appear minima that arise from balancing 
renormalizable versus non-renormalizable terms. A simple way to characterize them is to test 
the scaling v oc v^M in the large M limit, keeping all other Lagrangian parameters fixed (see 
|10j for more details). In contrast, for MSSM-like minima the VEV reaches a finite constant 
in this limit. We will use this criterion to characterize the type of vacua. 

4.4 Elect roweak constraints 

Lastly, to assess the viability of a given parameter point, we also estimate the contributions 
to the oblique parameters in order to select points that are in reasonable agreement with 



the EW precision constraints. As seen in Eq. (33), the heavy physics can induce tree- level 



contributions to the Peskin-Takeuchi T parameter [28] (see also Appendices |A] and |B| . If one 
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requires |T| ~ 0.2 and constrains Ci, C2 and C3 one at a time, one sees that for tan/3 > 1 
and M = 1 TeV, the strongest constraint is C2 < 0.1. Due to the and t^^ suppression 



in Eq. (33), Ci and C3 are less constrained for tan/? > 1. For instance, if tan/3 = 2 one has 
Ci < 1.3 and C3 < 0.15 (and are virtually unconstrained for tan (3 = 20). Notice also that there 
can be partial cancellations in Eq. ([33]). p| Thus, even if ci, C2 and C3 are of order one, provided 
M ~ 1 TeV, these contributions to T are not necessarily much larger than the experimental 
limit for a reference Higgs mass of ttlh^^i = 117 GeV. 

There are other contributions to T that can be comparable to the tree-level one -which is 
suppressed by v'^/M^- from loops of the MSSM Higgs fields (SM-like Higgs heavier than mn.^f, 
and mass splittings among the non-standard Higgses) as well as from custodially violating 
mass splittings in the (light) superpartner sector. The former give the following contributions 
to T and ^ [301 EI]: 

+ Apsmirrih) + ApgM(m^) - ApsmimH^J , (56) 
S""'''' ^ Y^{(^lL)'i^(m.,m^) + (r^f^j2F(mH,m^)-(r^|V^,)Mnm|± 

+ hZZ^\nml + HZZ^\nmjj-^-\nmj^^^^^ , (57) 
where a is the fine structure constant, i^ref is the reference SM Higgs mass used in the fit 



to the EW precision measurements, hVV/HVV with V = Z,W are defined in Eqs. (38) 



Vw±H^A = U'^w^H'fA + Vw^AH'f), with r]yy±H^A, Vw^AHT as defined in Eqs. (|40| [and similar 
definitions for the other ?7'^^'s], and 



Q 

[mf — 1712) (mf — \'^2 



^As shown in Appendix [a] Higgs triplets give ci,C2,C3 > 0. This corresponds to a positive tree- level 
contribution to T from triplets with zero hypercharge and negative from triplets with Y = ±1. Ref. [29\ 
points out that for lighter triplets it is also possible to have a cancellation in T arising from the interplay 
between the triplet mass and the /i-term. 
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where Hi = h, H. It should be noticed that the expressions (56) and (57) are only approximate 
since there can be logarithmically divergent terms in the loop computation of the T and S 
parameters, in the effective theory with higher-dimension operators, that scale like v'^/M^ 
(the quadratically divergent contributions to the gauge boson self-energies vanish by gauge 
invariance). However, due to the loop suppression, these contributions are expected to be 



much smaller than the tree-level one, given by Eq. (33), and are therefore negligible for points 
where the latter is within the experimental limits (the logarithmic enhancement is small for 
the scale of new physics we consider). 

As remarked in [32] when the SUSY particles are light there can be additional relevant 
contributions to the T parameter. These depend on parameters that do not affect directly the 
Higgs sector. As a result, we do not perform here a detailed fit to the EW data. However, we 
use Eqs. (56) and (57) plus the tree-level contribution, Eq. (33), to estimate whether a given 
point can be reasonably expected to agree with the precision constraints if appropriate SUSY 
contributions were added. In Ref. [32] it was found that such SUSY contributions are positive 



and easily as large as AT^^^^ ~ 0.2. Taking -0.2 < T*°* < 0.3 (95% C.L.)|^we therefore 
allow only points with —0.4 < T + T^'^^s < After this cut we find that all points in the 
scan satisfy —0.05 < S^''^^^ < 0.08, so that we do not impose any further cuts on S. 



4.5 Loop effects 

We have implemented our tree-level expressions for the spectrum and Higgs couplings in 
HDECAY v3.4 [3S]- This allows us to also take into account the QCD radiative corrections, 
that are known to be sizable. In addition, we include the radiative corrections derived from 
the 1-loop RG improved effective potential due to the supersymmetric particles, as computed 
in [37] . Loop contributions from the heavy physics that has been integrated out are suppressed 
by both a loop factor and by powers of M, hence they are expected to be negligible. 

As we will argue in the next section (see also Fig. [T]), already at tree-level the Higgs 
spectrum satisfies the bounds from LEP in large regions of parameter space. Hence, the 
motivation for taking rather large stop masses is absent in the extensions that we study. 
Keeping with the philosophy that the SUSY breaking scale is small compared to M, we 
consider a relatively light superpartner spectrum. For concreteness we take the superpartner 
soft parameters to have a common value Msusy = 300 GeV and = A;, = 0. |^ Thus, the 

^We use the code of Ref. [33J with mj = 173.1 ± 1.3 GeV/c^ [34;, Mw = 80.432 ± 0.039 GeV/c^ [ST, and 
mn^^t, = 117 GeV. 

^We denote the soft breaking masses by Msusy- We evaluate the scale inside the logarithms associated 
with SUSY loops at V^'^lt/sr + ~ 347 GeV. 
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SUSY loop contributions to the Higgs masses are modest, but the loop contributions to the 
Higgs couplings are more important and sensitive to the details of the SUSY spectrum |3Hl |39] . 
The point to remember is that the relevant loop-level effects can be fully computed given the 
MSSM superpartner spectrum, and are only mildly dependent on the details of the UV theory 
that gives rise to the EFT that we study. 



5 Numerical analysis: Results 



In this section we present the results for the Higgs spectrum and couplings that arise from a 



scan over the parameter region defined in Subsection |4.1[ We make sure that we concentrate 
on parameter points that are expected to be relatively insensitive to higher orders in the 1/M 
expansion, that they correspond to global minima, and that they can be in agreement with the 
EW precision constraints, given the uncertainties from the SUSY spectrum, as discussed in 



Subsections 4.2, 4.3 and 4.4 Since our choice of points depends on the "robustness" criterion 



described in Section 4.2, we start in Subsection 5.1 by commenting on the consequences of 
this prescription. We present these results in the tree-level approximation, so as to emphasize 
the effects of the higher-dimension operators. 

Next we present the results for some selected observables, such as the gluon fusion pro- 
duction cross section relevant at hadron colliders, and comment on certain "exotic" branching 
fractions. For these we include the supersymmetric radiative effects to the Higgs masses and 



couplings, as described in Subsection |4.5[ We will present a more complete study of the Higgs 
phenomenology in [ll], where we will also include the detailed bounds from LEP and the 
Tevatron on the Higgs spectrum. 



5.1 Sensitivity against higher-order corrections 



We described in Subsection 4.2| a simple prescription to estimate the sensitivity of a given point 
in parameter space against 0{1/M^) effects. In this subsection we illustrate the dependence 
on the allowed variation in v by requiring that the VEV change by less than 5%, 10% and 
20% if the next order corrections take a "typical" size, as estimated in Subsection 42_. As an 
example, we show in Fig. |2]m/i as a function of for tan (3 = 2 (left panel) and tan /? = 20 
(right panel). We have scanned over a total of 10^ points, but show only those points that 
correspond to a global minimum and that can be in reasonable agreement with the EW 
precision constraints, as explained in Subsection 4.4 For tan/3 = 2 (tan/3 = 20), we find 
that about 70% (80%) of the points in the scan correspond to global minima of the potential 
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Figure 2: Illustration of the sensitivity ofrrih against higher order effects in the 1/M expansion, 
as a function ofruA, for tan f3 = 2 (left panel) and for tan f3 = 20 (right panel) . The regions in 
(green, green + blue, green + blue + red) correspond to the requirement (6v/v < 0.05, 6v/v < 



0.1, Sv/v < 0.2), according to the prescription described in Subsection 4-2 The gray points 
are the additional points in the scan that do not obey any of these three requirements. The 
dashed blue line corresponds to the maximum tree-level value of rrih when the dimensionless 
coefficients of the higher- dimension operators are allowed to be as large as 1 (in absolute 
magnitude). The region of parameter space in the scan is described in the main text. 



y = Ken. + Kion-ren., defined by Eqs. and (12), while the rest are only local minima that 



we discard. The EW precision constraints further reduce the number of potentially viable 
points in the scan by 70% (80%). In Fig. |2]we exhibit how the number of points is further 
reduced by requiring 6v/v < 0.2 (green + blue + red regions), 6v/v < 0.1 (green + blue) 



and 5v/v < 0.05 (green region), following the prescription of Subsection 4.2 (in the figures 
we apply the requirements on 6tan(3 described in that subsection). Recall that this gives a 
measure of the sensitivity of a given point against higher orders in the 1/M expansion. The 
points shown in gray are the additional points that would change by 6v/v > 0.2 under the 



modification of Eq. (55), and without any restriction on 6tan(3 (i.e. in Fig. [2] we show all the 
points in the scan that are global minima and obey the EW precision constraints). 

We observe that the points with smaller mh are more easily affected by order correc- 
tions (for instance, they can lead to no EWSB after such a perturbation). This is not to say 
that there can be no models with small rrih, but only that points were rrih is heavier are rela- 
tively insensitive to those higher-order corrections. The reason is that larger rrih indicates that 
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both the 1/M and effects are contributing fully, without major accidental cancellations 



(see Subsection 4.2). The next order is then suppressed by order v/M, as expected (again, 
note that all our dimensionless couplings are at most one). Note also that for tan/5 = 20, 
the points with small itla are found to be rather sensitive to the higher-order corrections, and 
tend to be discarded using our prescription. Again, this is not to say that viable models with 
large tan/? and small do not exist, but only that their properties may not be correctly 
captured at the order we are working, so we choose not to concentrate on such cases. In 
the following, we restrict to points that satisfy 5v/v < 0.1, which should be interpreted as 
points for which the higher-order corrections introduce an uncertainty of at most order 10%. 
However, for many points the expected uncertainty should be smaller. 

In the figure we also show the maximal value of mh for the parameter region defined in 



Subsection |4.1[ This envelope was obtained by optimizing the values of the dimensionless 
model parameters so as to maximize rrih- The reason that the points in the scan itself do not 
reach such large values of mh is that there is a low probability that all the model parameters 
simultaneously attain the optimal values that maximize mh- 

5.2 Higgs Masses: Comparison to the MSSM 

We start by presenting our results for the Higgs spectrum, and comparing it to the MSSM 
one. We do this at tree-level only. Recall that the radiative corrections are "common" in 
the MSSM and in the effective theory under study, arising mainly from QCD and the MSSM 
superpartner sector. The observed differences can therefore be interpreted as arising directly 
from the heavy physics through the higher-dimension operators. We present in Figs. |3]j5] the 
results of the scan for m^-, mn and mH± as a function of m^, for tan/5 = 2 (left panels) and 
tan/5 = 20 (right panels). The tree-level MSSM curve is shown for the corresponding and 
tan (3. We also indicate which points correspond to sEWSB vacua (red points) and MSSM-like 



vacua (blue crosses), as described at the end of Subsection 4.5 We use the criterion described 
in [To], which is based on the fact that for large M with all other parameters fixed, the VEV 
in sEWSB vacua scales like v oc v^M. We see that for small tan 13 a large number of points 
are of the sEWSB type, while for large tan (3 most points correspond to MSSM-like vacua. 

As expected from our discussion of Section [3} the corrections to mh are more important 
at low tan/?. However, the scan shows that they can also be relevant at large tan/5. As 



already remarked in Subsection 5.1 there are significantly fewer points with small m/j, which 
is a consequence of the "robustness against higher-order corrections" criterion described in 
Section |4l 
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Figure 3: Lightest CP-even Higgs tree-level mass, rrih, as a function of itla, for tan (3 = 2 (left 
panel) and for tan (3 = 20 (right panel). The dashed blue line corresponds to the maximum 
value of mh when the dimensionless coefficients of the higher- dimension operators are allowed 
to he as large as 1 (in absolute magnitude). The region of parameter space in the scan is 
described in the main text. The solid green line corresponds to the tree-level MSSM result. 
Red points correspond to sEWSB vacua, while blue crosses correspond to MSSM-like vacua. 

It is also interesting that most points in Fig. [3] present significant deviations from tlie cor- 
responding MSSM values. Our parameter region includes the case where all higher-dimension 
operators vanish, and therefore includes the MSSM limit. However, in the scan it is unlikely 
that all of them become small simultaneously, which explains why there tends to be more 
points that exhibit important deviations in rrih with respect to the MSSM. In the large tan (3 
case, the overlap with the MSSM for sufficiently large rrih is possible mainly because many of 
the operators are l/tan/5 suppressed, and hence, at large tan (3 the number of relevant coeffi- 
cients that contribute to the departure of the Higgs spectrum from the MSSM one is smaller 
than at low tan [3. This implies that at large tan f3 there is a higher probability to effectively 
turn off all the effects from the higher-dimension operators, and reproduce the MSSM values 
for rrih- Our study simply encapsulates the picture of heavy physics of mass M (= 1 TeV) 
with couplings of order one to the MSSM Higgs sector. For instance, the light CP-even Higgs 
can become sufficiently heavy for the W^W~ and ZZ channels to be kinematically accessible, 
thus potentially allowing for a SM search of the light supersymmetric Higgs in the dilepton 
plus missing energy as well as the four-lepton channels. 



27 




Figure 4: Heavy CP-even Higgs tree-level mass, mn, as a function ofrriA, for tan j3 = 2 (left 
panel) and for tan/3 = 20 (right panel). The solid green line corresponds to the tree-level 
MSSM result. The region of parameter space in the scan is described in the main text. Red 
points correspond to sEWSB vacua, while blue crosses correspond to MSSM-like vacua. 




Figure 5: Charged Higgs tree-level mass, mH±, as a function of tua, for tan (3 = 2 (left panel) 
and for tan/? = 20 (right panel). The solid green line corresponds to the tree-level MSSM 
result. The region of parameter space in the scan is described in the main text. Red points 
correspond to sEWSB vacua, while blue crosses correspond to MSSM-like vacua. 
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Figs. [4] and [s] show the H and spectra as a function of the CP-odd Higgs mass for 
two values of tan (3, respectively. The deviations in tuh and mH± from the MSSM values are 
less dramatic in the large region, but they can be phenomenologically significant for lower 
values of rriA- Notice for instance that for ttla ~ 100 GeV, mn can be sufficiently large for 
the decay H AA to be kinematically open. Similarly, the decay W^A can be open. 

We will further comment on these decay channels in the next section. 



5.3 Higgs couplings to gauge bosons and fermions 



In Fig. |6] we present the tree-level couplings (normalized by the SM values) of the CP-even 
Higgses to Z pairs and up-type and down- type fermions, as a function of m^. The solid curves 
show the corresponding MSSM tree-level prediction, making it clear that large deviations from 
the MSSM can be induced via the higher-dimension operators. 

The plots in the upper row of Fig. [6] show the couplings of the CP-even Higgses to Z pairs 
for tan/5 = 2 (left panel) and tan/5 = 20 (right panel). We will refer to the Higgs state that 
has the largest coupling to the Z as "SM-like" (i.e. h is "SM-like" if \ghzz\ > IdHzzl while 
H is "SM-like" if \gHzz\ > \ghzz\)- Recall that the couplings to WW and ZZ are different 
as a result of the corrections given in Eq. (37). The difference appears at order and is 

of order a few percent when these couplings are sizable (e.g. for the SM-like Higgs). For a 
non-SM-like Higgs with a very small coupling to W^s and Z's, the relative difference between 
the two can be significant. We also note that in the present scenarios gf^zz + Qhzz need not 
add up to one, refiecting the fact that there are additional (heavy) degrees of freedom that 
have been integrated out. This effect also arises at order and we find that for the chosen 

parameters it can be as large as 7% (this can be seen more clearly in the large region of 
the figures). 

For large m^, h becomes SM-like to a very good approximation, as in the MSSM. However, 
we also see that in this limit H can have less suppressed couplings to Z pairs than in the 



MSSM. These are the properties anticipated in Eq. (53). At low wla and for tan/5 = 2, we 
see that H becomes SM-like to a better degree of approximation than in the MSSM, while h 
has more suppressed couplings to Z pairs than in the MSSM. There are also a large number 
of models with intermediate uia where both CP-even Higgses couple significantly to gauge 
boson pairs. For tan/5 = 20 the situation is similar in the large itla region. However, there 
are very few points at low itla (see Subsection 5.1 for an explanation). 

We also notice several features of the tih and hhh [ttH and bbH) couplings, which are 
relevant for gluon fusion induced processes. We consider first the case of low tan /5. The tth 
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Figure 6: Tree-level couplings of the CP-even Higgs states (normalized to the SM ones) to 
gauge bosons (upper plots), to top pairs (middle plots) and to bottom pairs (lower plots), 
as a function of ruA, for tan/3 = 2 (left panels) and taii/S = 20 (right panels). The blue 
crosses correspond to the couplings of h, and the red dots to those of H . The solid green 
lines correspond to the tree-level MSSM result. The region of parameter space in the scan is 
described in the main text. 
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Figure 7: Gluon fusion cross section at LO in Og, including light SUSY particle loops and 
normalized to the SM, for the light (upper panels) and heavy (lower panels) CP-even Higgses, 
and for tan/5 = 2 (left panels) and tan/? = 20 (right panels). The region of parameter space 
in the scan is described in the main text. We indicate by red crosses those points where h is 
SM-like (\ghzz\ > \gHzz\) and by blue dots those points where H is SM-like (\gHzz\ > \ghzz\)- 
The dashed line corresponds to the MSSM result. 
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coupling can present a small enhancement over both the SM and the MSSM of up to about 1.2 
when mA > 100 GeV (for smaller itla it can present an enhancement of up to 50% with respect 
to the MSSM). This would imply a factor of up to 1.4 with respect to the SM prediction for 
gluon fusion production relevant for searches at the Tevatron and LHC. In addition, given 
that there is no need for large radiative corrections to rrih from the squark sector to avoid the 
bounds on the light Higgs, if a light stop spectrum of order a few hundred GeV is present, 
additional relevant contributions to the gluon fusion process will be expected. These can lead 
to an additional enhancement of the reach in gluon fusion induced channels [IQ]. The bbh 
coupling is suppressed both with respect to the SM and the MSSM when rriA > 200 GeV, and 
reaches a value at most equal to tan (3 in the low rriA region (as in the MSSM). However, notice 
that there are a number of models with suppressed couplings to bb pairs for rriA < 140 GeV. 

For low tan/5 the contribution to gluon fusion is governed by the top Yukawa coupling, 
and therefore there is a net enhancement with respect to the SM for many models, as seen in 
the upper left panel of Fig. [7| where we present the gluon fusion cross section for h production 
normalized to the SM one, as a function of nih- In this figure we include the radiative effects 



to the Higgs masses and couplings assuming a SUSY spectrum as described in Subsection 4.5 



Most of the models with m/j < 115 GeV are excluded by LEP, while the Tevatron excludes 
some of the models with rrih in a window around 170 GeV. A detailed analysis of all the 
cross sections and branching fractions to determine the allowed models will be presented 
elsewhere [Hj. In Fig. [7]we indicate by red crosses the models where h is "SM-like", as per 
the definition at the beginning of this subsection, and by blue dots those models where H is 
"SM-like". 

The tiH coupling at low tan/3, shown in Fig. [6| is found to be generically suppressed 
with respect to both the SM and the MSSM, except for the region rriA < 200 GeV where 
some enhancement may be possible. The bbH coupling is generically enhanced with respect 
to the SM and even with respect to the tan/? enhanced value in the MSSM for large rriA- 
However, the bottom loop contribution to the gluon fusion process is subdominant, as seen in 
the lower left panel of Fig. [7| where a net suppression with respect to the SM is observed for 
niH > 200 GeV, governed by the suppression in the tiH coupling. In this figure we also see 
an enhancement with respect to the SM for itlh < 200 GeV, which reflects the enhancement 
of the ttH coupling for < 200 GeV mentioned above, in addition to the enhancement due 
to light stop contributions. 

^We compute (j{gg h) /a^^{gg h) ^ T{h — > gg) /T^^{h gg), which holds at leading order 
in Us [SJ 1121 The values of T{h gg) and T^^{h gg) are calculated at LO using a version of 

HDECAY f35^ , modified to include the tree-level expressions in the presence of the higher-dimension operators. 
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Figure 8: Tree-level couplings of the CP-odd Higgs to fermion pairs as a function of niA, 
for tan/5 = 2 (left panel) and tan/3 = 20 (right panel), and normalized to the SM value, 
gruf /2mw . The blue crosses correspond to the couplings of A to ti, and the red dots to those 
of A to bb. The solid green lines correspond to the tree-level MSSM result. The region of 
parameter space in the scan is described in the main text. 

At large tan f3, there is no significant variation in tth with respect to the SM or the MSSM 
for niA > 200 GeV, but there is a small suppression for 100 GeV < niA < 200 GeV, and there 
can be an enhancement with respect to the MSSM for rriA < 100 GeV. The bbh coupling is 
enhanced with respect to the SM and the MSSM in a large number of models, and achieves 
the largest values for rriA < 100 GeV, although it is smaller than the tan j3 value that occurs 
in the MSSM. For mA > 200 GeV there are many models where the bbh coupling is strongly 
suppressed. As shown in the upper right panel of Fig. [7], in many models the regions of 
enhanced bbh coupling lead to a relevant enhancement of the gluon fusion cross section with 
respect to the SM one. 

At large tan (3, the tiH coupling can have an enhancement with respect to the MSSM in 
the region 100 GeV < < 200 GeV, but is generically suppressed with respect to the SM. 
The bbH coupling presents the familiar tan /3 enhancement of the MSSM at large ttla, although 
there could be some suppression in the intermediate region 100 GeV < < 200 GeV. At 
low rriA we see a few models with a significant enhancement of the bbH coupling with respect 
to the SM and the MSSM. 

In the large tan (3 region both top and bottom loops contribute to the gluon fusion cross 
section with different weights depending on the ratio of the fermion to Higgs masses and the 
specific (large) value of tan/3. The lower right panel of Fig. [T] shows that for values oi mn 
larger than 200 GeV there is always a suppression with respect to the SM gluon fusion cross 
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Figure 9: Branching fractions for H AA (left panel) and H AZ (right panel). The 
dashed line corresponds to the MSSM result. SUSY and QCD radiative corrections are in- 
cluded. The plots are for tan/3 = 2, M = 1 TeV, /i = m, = 200 GeV, Msusy = 300 GeV, 
At = Ah = 0, and a scan over the ranges defined in Section We indicate by red crosses 



those points where h is SM-like (\ghzz\ > \gHzz\) md by blue dots those points where H is 
SM-like (\gHzz\ > \9hzz\)- 

section for H, similar to the MSSM case. 

The couphngs of the CP-odd and charged Higgs bosons differ from those of the MSSM 



only at order due to the corrections to their kinetic terms [see comments after Eq. (39)]. 

These deviations are far less significant than for the CP-even Higgs states. As an example, 
we show in Fig. [8] the couplings of the CP-odd Higgs to up-type and down-type fermion 
pairs, which follows closely the tan/? enhancement /suppression familiar in the MSSM [see 



also Eq. (35) 



We postpone a detailed discussion of all the relevant branching fractions and the conse- 
quences for Higgs searches at the Tevatron and the LHC to Ref. [H]. Here we comment on 
a number of selected "exotic" channels that can motivate new search strategies, in particular 
regarding the heavy CP-even and charged Higgs bosons. These are mostly related to the 
distortion of the Higgs spectrum with respect to the MSSM which can lead to the opening 



of new decay channels, as mentioned at the end of Subsection 5.2 For example, we show in 
the left panel of Fig. |9]the branching fraction for H —>■ AA, which can be significant for mn 
up to about 200 GeV. It can be in fact much larger than in the MSSM, where such values of 
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Figure 10: Branching fractions for H —>■ hh for tan (3 = 2 (left panel) and tan/3 = 20 (right 
panel). The dashed line corresponds to the MS SM result. SUSY and QCD radiative corrections 
are included. We indicate by red crosses those points where h is SM-like (\ghzz\ > \gHzz\) 
and by blue dots those points where H is SM-like (\gHzz\ > \9hzz\)- 

mil are already near the decoupling limit so that such decays are highly suppressed by phase 
space. Similarly, there can be important branching fractions for H AZ, as seen in the 
right panel of Fig. |9l Note that for most of these points it is H that is SM-like. 



We also show in the left panel of Fig. 10 that for low tan (3 there is a generic suppression 
of the H hh channel compared to the MSSM, except above the tt threshold where some 
enhancement is possible (the pronounced dip at mn ~ 108 GeV in the MSSM curve is due 
to an accidental cancellation of the Hhh coupling; in this small window H decays mostly into 
6's). The above suppression can be very significant for 170 GeV < mii < 250 GeV, reflecting 
a relatively heavy h so that the channel is closed. At large tan (3, however, there are large 
regions where BR(i7 hh) is enhanced with respect to the MSSM, and in fact H ^ hh can 



we 



be a significant decay channel in this case (see right panel of Fig. 10). 

Finally, we comment on certain charged Higgs decays. In the upper row plots of Fig. 11 
present the branching fractions for W^h {W^*h)., which show a significant suppression 

with respect to the MSSM at low tan (3 (except above the ti threshold, where some enhance- 
ment is possible), while for large tan/5 there is an enhancement compared to the MSSM in a 



large number of models. In the lower row plots of Fig. 11, we see that the branching fraction 



for W^A {W^*A), can be sizable at low tan/3 and for 100 GeV < m^i < 180 GeV, 
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Figure 11: Branching fractions for — > W^^*^h (upper panels) and — > (lower 
panels), /or tan /5 = 2 (left panels) and tan f3 = 20 (right panels). The dashed line corresponds 
to the MSSM result. SUSY and QCD radiative corrections are included. We indicate by red 
crosses those points where h is SM-like (\ghzz\ > IohzzI) o-nd by blue dots those points where 
His SM-like (\gHzz\ > \9hzz\)- 
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while it remains small for large tan jS (although it can still be enhanced compared to the 
MSSM). Since at low tan/?, BR(i7^ — ^ 14^^ A) can be close to unity, one may therefore ex- 
pect to produce a large number of CP-odd Higgs bosons in top decays. Note also that a 
large fraction of the points with such a property present the "inverted hierarchy" where H is 
SM-like. 

6 Conclusions 

We have considered a large class of supersymmetric scenarios with physics beyond the MSSM 
that couples appreciably to the MSSM Higgs sector. Our main assumption is that the degrees 
of freedom beyond the MSSM are heavier than the weak scale, and that their SUSY mass split- 
tings can be treated as a perturbation. We call this approximately supersymmetric threshold 
M. As a result, a model-independent analysis can be set up, based on an approximately 
supersymmetric effective theory of the MSSM that includes higher-dimension operators sup- 
pressed by powers of 1/M. These higher- dimension operators can encapsulate different types 
of physics beyond the MSSM, such as singlet or triplet Higgses, heavy Z' and W, etc. (we 
illustrate the detailed connection in Appendix [A|). 

We argued, based on the structure of the induced Higgs quartic couplings, that both the 
leading and next-to-leading order in the 1/M expansion can be phenomeno logically relevant, 
and computed the Higgs spectrum and couplings up to this order. This included a careful 
treatment of degenerate cases and the inclusion of kinetic term renormalization that contains 
information about the mixing of the light and heavy degrees of freedom. The most important 
effects of the new physics enter through the angle a that characterizes the mixing in the CP- 
even sector, but we have systematically included all the effects to order This allows us 
to include in our analysis the recently discussed sEWSB vacua ^0\, which depend crucially 
on certain dimension-6 operators. 

We were especially careful to single out points in the effective theory that can be expected 
to be reliably described by the EFT. We also made sure that these points correspond to global 
minima of the effective potential (we did not consider the possibility of long-lived metastable 
minima). In addition, we took into account the EW precision constraints, making sure that 
the study points can be in agreement with precision tests when possible effects from squarks 
and sleptons not directly related to the Higgs sector are included. 

The large class of SUSY models presented in this study has in part already been explored 
by various Higgs searches at LEP and the Tevatron. We will present these constraints in the 
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accompanying paper yL4j. Similarly, we defer a more complete study of the Higgs collider 
phenomenology to that work. Here we simply pointed out a few interesting features that 
arise from our study: a generic enhancement of the gluon fusion production cross section of 
the SM-like Higgs, the presence of new channels with significant branching fractions, such as 
H A A and — > W^A, and the possible suppression of decay modes such as H ^ hh. 

We find it interesting that a study of the light Higgs sector can indirectly reveal the 
presence of new physics that either may not be within reach, or may not be easy to produce, 
as was also emphasized recently in |l3l 01]. The measurement of the Higgs spectrum and 
observation of some of their decay modes, together with the observation of some of the SM 
superpartners, can give a striking evidence for a more complicated structure beyond the 
MSSM. 
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A UV completions 

Here we consider possible UV completions to exemplify how the operators of the effective 
theory can be generated. We consider the addition of singlets, of vector-like triplets {Y = ±1), 
of a triplet with a Majorana mass and Y = 0, and an SU{2) gauge extension under which 
Hu and Hd transform as doublets. This illustrates how the higher-dimension operators on 
which the EFT analysis presented in the main text is based can arise, and also shows that 
the coefficients of these operators can be fairly uncorrelated if the extension is sufficiently 
complicated. This justifies our approach of scanning over these coefficients without taking 
into account the correlations that may arise in a given UV model. 
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A.l MSSM + Singlet 

The simplest model we consider shares structural features with the NMSSM. In the NMSSM, 
a new singlet S is added to the MSSM to solve the /i problem, through a VEV for S and the 
superpotential interaction XsSHuH^. In this model, the low energy theory contains additional 
scalar states from S, and after EWSB the singlet state mixes with the Higgs. However, there 
is a different region of parameter space, where the singlet S has a mass somewhat larger than 
the electroweak scale, and is integrated out of the effective theory. Although the fi problem is 
not addressed in this limit, the low energy theory contains modifications in the Higgs sector 
that can help solve the little SUSY hierarchy problem. 
The superpotential and Kahler potential are 

W = fiH^H,+ ^MsS^ + XsSH^H,- X (^a.fiH^H. + ^a^MsS^ + asXsSH^H,^ ,(58) 

K = Hle''H^ + HyHd + S^S-X^x{h^HlHd + h2HiHu + hS^S) , (59) 

where the and hi are dimensionless, and an e tensor is understood in the contraction 
HuHd = Hu^Hd- One can consider a cubic coupling kS"^, but it does not contribute to the 
effective theory up to the order we have analyzed, hence we do not write it explicitly 



in Eq. (58). The parameters 6i, 62 and ai map into fnH^-, n^Hu ^^e Bfx term, respectively. 



Integrating out S at tree- level induces nonzero cui, ai, C4, and 74 in Eqs. ((T])-([7]) as follows: 

M = Ms , ui = -Xl , ai = a2 - 2a3 , 

C4=\Xs\ , 74 = 02-03, P4= 102-03] -63, 

while all other EFT coefficients vanish. 

A.2 MSSM + SU{2)l Higgs Triplets 

Consider an extension with two SU{2)l triplets, T and T, with hypercharges Y = —1 and 
Y = +1, respectively. The superpotential and Kahler potential are: 

W = fiH^H, + MTTT+^XTH^TH^ + ^XfHafHa 

- X (^iiiH^Ha + a^MrTf + ^asXrH^TH^ + ^a.XfHdTHa^ , (61) 

K = HyH^ + Hy^H^ + T^e^^T + T^e^^T 

- XXt (biH^Hd + b2HlH^ + hT^T + b^f^f'^ . (62) 
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M = 






— -At At , 




= a2 - 


- O3 — O4 , 


Cl = 




71 


= 02 — 04 , 




= ^2 


- 04!^ - 63 


C2 = 


^IAtI , 


72 


= 02 — 03 , 


/32 


= ^2 


- O3P - 64 



An epsilon tensor is understood in the contractions HyTHu = HueTHu, etc. The parameters 
bi, 62 and oi map into rn'jj^, and the 6/i term, respectively. Integrating out the triplets, 
one can write the effective Lagrangian in terms of the operators defined in Eqs. ([T])-([7]), 
withF] 



(63) 



and all other dimensionless coefficients vanishing. 

If on the other hand one considers an extension with a single SU{2)i triplet with hyper- 
charge y = and a Majorana mass Mt, the superpotential and Kahler potential are: 

W = fiHuHd+^MrT' + XTH^THd-X (^aifiH^Hd + ^a2MTT' + a3\THuTHd^ (64) 

K = Hle^''Hd + Hle^^Hu + T^e^^T-XX^(b^H\Hd + b2HlHu + b:,T^T^ , (65) 
and after integrating out the triplet, one obtains the effective theory of Eqs. ([T|-([7]), with 



(66) 



while all other EFT coefficients vanish. 
A. 3 SU{2) Extensions 

The higher-dimension operators induced by integrating out a massive f/(l)' were worked out 
in Ref. [Bj. Here we consider a product gauge group SU {2)i x SU {2)2 x U{1)y with gauge 
couplings (71, g2 and g' , respectively. We concentrate on the SU {2)i x SU {2)2 factor since the 
U{1)y factor enters in a trivial way. The MSSM Higgs superfields are assumed to transform 
like (2,0) under SU {2)i x SU {2)2. The gauge group is broken down to the diagonal by a S 
field transforming like (2, 2) under SU {2)i x SU{2)2, which gets a nonzero VEV, (S) oc cr^. In 

lOHere we use the identity / d'^OA^e^^A =\jd^O (We^^H)^, where A = A^X" with A" = Her^H, and 
X"' are the SU{2) generators in the adjoint representation while H is in the fundamental representation of 
SU{2). One also has Jd'^e A'' e^^ A ^ Jd^O {^{Hle^^ Hu){Hle^^ Ha) - i|iJ„ei?rf|2} for A" = Huer^Hd. 



M = 


Mt , 




- -4AT , 


ai = 


= 02- 


203 , 




cs = 


^lA P 


73 


= 02 — 03 , 




= 02 


- 03! 


-63 


C4 = 


^A P 


74 


= 02 — 03 , 




= 0-2 


|2 

- asl 


-63 
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this case, the kinetic term for S, TT{e'^92V^J:^e'^3^^^J:), leads to the structure TT{e'^32V2^2giVi-^^ 
which contains a mass term for the hnear combination V = {giVi + (72 V2)/ gf + 92- The 
orthogonal combination, V = {g2Vi — giV2) 



+ gi remains massless and can be identified 
with the SM W gauge bosons. Therefore, we take as the starting point the Kahler potential 



Ml, 



,2g2V2^2giVi 



(67) 



where g^^^ = gV + gV, g^V^ = {gl/glW - gV and g = gm/Vgf + gl is the SM SU{2)l 
gauge coupling while g = gf/ ^/ gj + gi is the coupling of the massive V. The term of second 



order in V in the last term of Eq. (67) identifies My as the gauge boson mass. In order 
to integrate out V to order 1/My,, it is sufficient to keep terms up to quadratic order in V 
in the "mass term", while keeping up to linear order in the terms not proportional to My,. 
We can also assume that V is in the Wess-Zumino gauge so that we can expand to quadratic 
order in V everywhere. Note that the term proportional to My, contains terms in addition 
to the pure mass term for V: 



M^, 



2igl + gl 



^2[(gi/gf)gV'-gV]^2{gV'+gV) 
1 



3_ 

3 



(y-V'^f _ (y^V) (K'V^)j \ . (68) 



The terms quartic in the vector multiplets are essential for maintaining the low-energy gauge 
invariancep] The equation of motion for V , to order 1/My,, gives 

.2 X 1 



V 



la 



^ i=u,d 



6 



ybyb j ^ _gya (^3 ^ 2^ VV^ + 2gV\^) \ H, 



(69) 



where there is no summation over a (but there is over 6), and are the SU{2) genera- 
tors in the fundamental representation. Replacing back in the action leads to the effective 
Kahler potential 



K, 



eff 



(TO) 



where we restored the terms in V into the exponential form, and the hypercharge vector 
multiplet can be put back in a trivial way. Therefore, the massive W induces the operators 
of Eqs. ([l|-(l7|) with coefficients 



Cl 



C2 



C3 



C4 = 2g' 



(71) 



-'^-'^ These terms are not present in the U{iy case. For instance gauging only the U{1)2 C SU {2)2 generated 
by r^, i.e. setting V2 = V2 = 0, the last two terms in Eq (68 1 cancel out. 
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M = My/, and all other EFT coefficients vanishing. 

There are also excitations of the S field (a SM singlet and triplet) with masses proportional 
to (S). These may be split by SUSY breaking and generate further contributions to the 
dimension-6 operators in the effective Kahler potential. These were considered in Ref. |1] in 
the opposite limit that interests us here, namely when SUSY breaking is comparable to the 
SUSY preserving contribution. 

For completeness, we also record here the SUSY-preserving operators induced by integrat- 
ing out a massive U{iy gauge field, i.e. integrating out V from 

where g' is the f/(l)' gauge coupling and Qu,d are the U{iy charges of Hu,d, respectively (for 
simplicity, we have omitted the SM gauge factors, which are implicitly understood). Unlike 
in the non-abelian case, the C4 operator is not induced, while 

ci = -AQlg[' , C2 = -4Qlg'^ , C3 = -AQ^Qdg" • (73) 

From the examples above, it should be clear that MSSM extensions including singlet and 
triplet Higgses, plus abelian or non-abehan gauge factors can generate all of the operators 
we considered in the main text, with essentially arbitrary coefficients, except for the Cq and 
C7 operators. We have nevertheless included the latter operators in our phenomenological 
analysis, since they are allowed by supersymmetry. 

B Comments on Custodial Symmetry 

Neglecting the Yukawa couplings, the Higgs sector in the MSSM can be written in terms of a 
chiral superfield 



Ha Hi 



as 



W D ^iHueHd = ^Tre^$^e$ , 

where W and B are the SU{2)l and U{1)y vector superfields, while e is the antisymmetric 
2-index tensor. This shows that in the limit g' ^ (and neglecting Yukawa couplings), the 
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theory is invariant under SU{2)l x SU{2)ii where $ — > Ul^U\^. The SU{2)l factor is gauged 
and can be complexified. The global SU(2)ji implies that in the limit Vu = Vd the theory has 
a custodial SU{2)l^^ global symmetry. This symmetry has a somewhat limited use since it 
is broken away from tan/5 = 1. 

However, at tree-level (and for g' = 0) the Higgs (scalar) sector of the MSSM exhibits an 
SU{2y^'^^^ X SU{2)^°^'^^ symmetry such that the Higgs scalar components 



transform as (j)u,d UL4>u,dUR- In fact, if the Higgsinos are taken to be singlets under the 
SU(2)ji transformation, only the Higgs-Higgsino-gaugino couplings (apart from the superpo- 
tential Yukawa terms) violate this global symmetry. Thus, for arbitrary expectation values, 
Vu and Vd, contributions to the T parameter enter only at one loop level. 

The fact that in the MSSM the tree-level gauge boson masses satisfy the SM relation p = 1 
can be seen directly from the identities Hi {2gWfHu = Tr (2^iy)2$„ and {2gWfHd = 
TT^li{2gWy^d, which lead to the gauge boson mass terms in the Kahler potential af- 
ter replacing the chiral superfields by their scalar components, ^u,d <Pu,d- This obser- 
vation also implies that higher-dimension operators such as [Hle^^ Hu)(HueHd + h.c.) or 
[H^e^^ Hd){Hu€Hd + h.c), i.e. those in Eq. (5), lead to p = 1 at tree-level. On the other 
hand, the higher- dimension operators of Eq. (4) contribute to the gauge boson mass terms 
through the linear term in e^^ = 1 + 2V + ■ ■ ■ , and these terms do not respect a custodial 
symmetry. This was seen explicitly in Eq. (33). As shown above, an exception is the operator 
{Hi e'^^Hu + H\ e^^HdY, in the hmit v^ = Vd- This operators is naturally induced by massive 
W primes, as shown in |A.3 



C Charginos and Neutralinos 

Although our focus is on the scalar Higgs sector, neutralinos and charginos can provide im- 
portant decay channels, and the effects of the higher- dimension operators on their masses 
play an important role. We have implemented in HDECAY the corresponding mass formulas 
to order 1/M, as computed in [10]. At order there are additional contributions to 

the mass matrix as well as to the kinetic terms, that are needed to compute the physical 
chargino/neutralino masses. Given that these states tend to be near the experimental limits, 
it would be interesting to compute these next order corrections, but we postpone it to future 
work. For completeness, we collect the mass matrices to order 1/M. 
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The chargino mass matrix is 



M2 V2- 



V2mwSi3 yu (1 - PS2/3) J \ 



w- 



(74) 



where p = uJiv"^ / (4/iM) takes into account the effects from the heavy physics. The neutrahno 
mass matrix is 



nizSwCfs mzswsp \ ( B \ 

1^3 



/ Ml 

M2 mzcwcp -nizCwSfs 

-mzSwC(3 mzCwCf3 '^f^ps} -p (1 - 2ps2 

\ nizswsp -nizcwsp -p (1 - 2ps2p) '^ppc} J 



(75) 



where cy/ stands for the weak-mixing angle cosOy/- Mi and M2 are the SU{2)l and f/(l)y 
gaugino soft breaking parameters. 



D Higgs couplings 



The couphngs between the Z gauge boson and two Higgs fields, defined in Eq. (40), are given 
by 

VZhA = Ca-f3 1 1 ~ ^ [(^1 + ^1) ^7 ~ ^I'^t] ^ ^\ ^^^^ ^ ~ "^"-^ 

+ c^Sr]zhA - s^6t]zha , 



+ c^Sr]zAh - s^6r]zAH 



(76) 



1 



VZHA = Sa-f3 I ^ ~ 2 ^^^-^ ^-^^ + ^1'^'^] j 2 ^^^^ + Ei)s^ + BiC^] Ca-f3 

+ c^Sr]zHA + SjdrjzhA , 

VZAH = Sa-f3 I ^ ~ ^ 1(^1 + ^7 + ^1-57] | ~ ^ [(^1 + ^1)^7 + BiC^] C^-IS 

+ c^5rjzAH + s^5rjzAh , 
= 1 - Fi + 5r]zH+H- , 



ZH+H- 
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where 7 is defined by Eq. (29), Ai, Bi, Di, Fi and Ei are defined in Eq. (27), and 



SVZhA = 


IVl 


L 


Sr]zAh = 


1 VI 


^/3^4 ( 

L 


5r]zHA = 






SrjZAH = 




^/3^4 ( 


St]zh+h- = 







(77) 



C3 (1 + 3s/ - 6t/ + tp') + tp' 



C6 {Cl + C2)tp +C7tp 

C2W 



Similarly, the couplings between the W and two Higgs fields, defined in Eq. (40), are given 
by 



Ca-/3 1 1 - ^ [(^1 + ^l) ^7 - ^1^7] 



1 



[{Di + Fi) - BiC^] Sa- 



■/3 



7J '^a—P 



+ C^5'qw±HTh - S^5riw±HTH 



(78) 



[(Ai + Fi) + BiC^] C„_ 



VW^H^A = 1 - - (Fi + Fi) + 5?7h/±//TA , 



Vw^AH^ = 1 - - (Fi + Fi) + 5rjw±AHT 
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where 



Sr]W±H^h = ^JpSl{c3{Sa + tp^Ca) +tp^ {citp^Sa + C2Ca)} , 

6vw±HH^ = ^sli^^C3[{3 + s/-9t/)c^+{t-\9-s/)-3t/)s^] (79) 

+ 3tp^ (cit^^Ca - C2Sa) | , 

^2 

5W±AHT = ^4{4^/(ci + C2) + C3(l + 3s^'^^-6t/ + t/)} . 
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